Issues with J-dependence in the LSDA+£7 method for non-collinear magnets 
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We re-examine the commonly used density functional theory plus Hubbard U (DFT+U) method 
for the case of non-collinear magnets. While many studies neglect to explicitly include the exchange 
correction parameter J, or consider its exact value to be unimportant, here we show that in the case 
of non-collinear magnetism calculations the J parameter can strongly affect the magnetic ground 
state. We illustrate the strong J-dependence of magnetic canting and magnetocrystalline anisotropy 
by calculating trends in the magnetic lithium orthophosphate family LiMPC^ (M = Fe and Ni) and 
difluorite family MF 2 (M = Mn, Fe, Co and Ni). Our results can be readily understood by expanding 
the usual DFT+17 equations within the spinor scheme, in which the J parameter acts directly on 
the off-diagonal components which determine the spin canting. 
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Density functional theory (DFT) within the local den- 
sity (LDA) and generalized gradient (GGA) approxima- 
tions is widely used to describe a large variety of materi- 
als with good accuracy. The LDA and GGA functionals 
often fail, however, to correctly reproduce the properties 
of strongly correlated materials containing d and / elec- 
trons. The LDA+£7 approach - in which a Hubbard U re- 
pulsion term is added to the LDA functional for selected 
orbitals - was introduced in response to this problem, and 
often improves drastically over the LDA or GGA. Indeed, 
it provides a good description of the electronic properties 
of a range of exotic magnetic materials, such as the Mott 
insulator KCUF3 1 and the metallic oxide LaNi02 2 . 

Two main LDA+U schemes are in widespread use 
today: The Dudarev 3 approach in which an isotropic 
screened on-site Coulomb interaction U e ff — U — J is 
added, and the Liechtenstein 1 approach in which the U 
and exchange (J) parameters are treated separately. The 
Dudarev approach is equivalent to the Liechtenstein ap- 
proach with J = 4 . Both the effect of the choice of 
LDA+£7 scheme on the orbital occupation and subse- 
quent properties 5 , as well as the dependence of the mag- 
netic properties on the value of U 6 , have recently been 
analyzed. There has been no previous systematic study, 
however, of the effect of the J parameter of the Liechten- 
stein approach in non-collinear magnetic materials. Here 
we show that neither the approach of not explicitly con- 
sidering the J parameter (as in the Dudarev implemen- 
tation), nor the assumption that its importance is bor- 
derline - a common approximation is to use J ~ 10% U 
without careful testing - within the Liechtenstein imple- 
mentation are justified in the case of non-collinear mag- 
nets. We demonstrate that in the case of non-collinear 
antiferromagnets, the choice of J can strongly change the 
amplitude of the spin canting angle (LiNiPC^) or even 
modify the easy axis of the system (LiFePC^ and FeF2), 
with consequent drastic effects on the magnetic suscep- 
tibilities and magnetoelectric responses. 

First we remind the reader how the U and J parame- 
ters appear in the usual collinear spin LSDA+/7 formal- 



ism. The LSDA+Z7 reformulation of the LSDA Hamilto- 
nian is usually written as: 

HlSDA+U = H LSD A + Hu , (1) 

whith 

mi ,m,2 

where P is the projection operator, a is the spin index, 
and (on a given atomic site): 

T/t(l) = 

7712 ,77li 
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+ (^3,2,4 - ^*1,2) 4™ + J)S lt2 (3) 

Here ^3,2,4 = (™>urn 3 |V^ e 1)m3>m2>m4 1 ra 2 , m 4 > are the 
elements of the screened Coulomb interaction (which can 
be viewed as the sum of Hartree (direct) contributions 
Vfl 2 4 and Fock (exchange) contributions V£% 4 2 and 
n?j are the d-orbital occupancies. 

In the case of non-collinear magnetism, the formal- 
ism is extended and the density is expressed in a two- 
component spinor formulation: 

P =(P]1 P] 1 ,) = ( n + m * m x -im y \ 

\p^ 1 p++ J yn x + trriy n — m z J v 7 

where n is the charge density and m a the magnetiza- 
tion density along the a direction (a = x,y,z). Using 
the double-counting proposed by Bultmark et a/. 7 , the 
LSDA+/7 potential is then also expressed in the two- 
component spin space as: 

where and are equal to Eqs.3 and 

« = £ (-KIa,2 + «& at) (6) 
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For collinear magnets, only and (Eqs. 3) are 
relevant since and are equal to zero, and J affects 
the potential mainly through an effective U—J. However, 
in the case of non-collinear magnetism, the and 
and hence the and (Eqs. 6) are non-zero. Then 
it is clear from Eqs. 6 that J acts explicitly on the off- 
diagonal potential components. 

Next, we show the effect of the choice of J parameter 
in the family of lithium orthophosphates, LiMP04 (M 
= Ni and Fe) and in the family of difluorites MF 2 (M 
= Mn, Co, Fe and Ni). The orthophosphates crystallize 
in the orthorhombic Pnma space group with C-type an- 
tiferromagnetic (AFM) order. The difluorites crystalize 
in the tetragonal P^/vnnm rutile structure with AFM 
order. We performed calculations within the Liechten- 
stein approach of the DFT+Z7 as implemented in the 
VASP code 8,910 with U and J corrections applied to the 
3d orbitals of the M cations. In all cases we relaxed the 
atomic positions until the residual forces on each atom 
were lower than 10 /ieV/A at the experimental volume 
and cell shape reported in Tab. I, taking into account 
the spin-orbit interaction. We found good convergence 
of the non-collinear spin ground state with a cutoff en- 
ergy of 500 eV on the plane wave expansion and a k-point 
grid of 2 x 4 x 4 for the orthophophates and 4 x 4 x 6 for 
the difluorites. 




(b) 


500 




400 


> 


300 


CD 


200 




100 




i— 
CD 





C 
LU 


-100 




-200 




-300 



(c) 
> 

CD 

LU 
< 
O 



U=5, J=0 — e- 
U=5,J=1 
U eff = 4 
U=5, J=1 (modified) — eh 




2 4 6 8 10 

canting angle (°) 




a 



Ref. 



LiFeP0 4 10.332 6.010 4.692 11 

LiNiP0 4 10.032 5.854 4.677 12 

NiF 2 4.650 4.650 3.084 13 

FeF 2 4.700 4.700 3.310 14 

MnF 2 4.650 4.650 3.084 15 

CoF 2 4.695 4.695 3.179 16 

TABLE I. Experimental cell parameters (A) used in the sim- 
ulations of LiMPCU phosphates and MF 2 difluorites. 



FIG. 1. (a) Calculated LSDA+[7 canting angle of LiNiP0 4 
versus J for U = 5 eV. The experimental value of the canting 
angle is equal to 7.8° 17 . (b) Energy versus canting angle in 
LiNiP0 4 for U = 5 eV and J = eV (red circles), U = 5 eV 
and J — 1 eV (blue triangles), U e ff = 4 eV (green crosses) and 
U — 5 eV and J — 1 eV but by fixing J — eV in Eqs. 6 (pink 
squares). The zero energy reference is chosen at zero canting 
angle, (c) Magnetocrystaline anisotropy energy (MCAE) be- 
tween the a and b orientations of the magnetic moments of 
LiFeP0 4 . The experimental b orientations is taken as energy 
reference. 



First, we focus on LiNiPCU, which is known experi- 
mentally to be C-type AFM, with an easy- axis along the 
c direction and a small A-type AFM canting of the spins 
along the a direction (C Z A X ground state with mm , m 
magnetic point group) 17 . Performing calculations within 
the LSDA+£7 method with J = 0, we find that we cor- 
rectly reproduce the C Z A X ground state with a rather 
small U sensitivity of the magnetocrystalline anisotropy 
energy (MCAE) and the spin canting; this finding is 
consistent with a previous report using the GGA func- 
tional 18 . However, our calculated canting angle of 1.6° 
for U = 5 eV and J = eV severely underestimates 
the experimental value of 7.8° 17 . In Fig.l (a) we show 
the evolution of the canting angle with J at U = 5 eV. 
We find that the canting angle is extremely sensitive to 
the value of J - in fact it is cx J 3 - changing from 1.6° 
at J = eV to 7.8° at J = 1.7 eV. To reproduce the 
experimental value of the canting angle we need to use 
the rather large J value of 1.7 eV. The dependence of 
the canting angle on J is consistent with Eqs. 6, as the 



off-diagonal elements and are non-zero when the 
spins cant away from the easy axis. 

In Fig. 1 (b) we report the energy versus the canting 
angle in LiNiPC^ for U = 5 eV and different values of J. 
We see that as J is increased from J = eV to J = 1 eV 
(red circles and blue triangles) the minimum of the en- 
ergy shifts to larger canting angle, with a stronger gain 
of energy with respect to the uncanted reference. When 
performing the same calculation with U e f / =4 eV (green 
crosses in Fig. 1) we obtain results that are very similar 
to the case U = 5 eV and J = eV, which is formally 
equivalent to the Dudarev approach with J7 e // = 5 eV. 
These comparisons confirm that varying U has a minimal 
effect on the canting angle in LiNiPC^ and also that the 
use of the Liechtenstein treatment of J is extremly impor- 
tant. To further confirm the direct relationship between 
the spin canting and the J parameter, we performed the 
same calculations with U = 5 eV and J = 1 eV but we 
artificially fixed J = eV only in Eqs. 6 (pink squares in 
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Fig.l (b)). We clearly see that the energy versus canting 
angle is strongly affected by this modification and in fact 
the canting is almost removed. 

Similar J dependence of the canting angle was also re- 
ported previously for Ni 2+ in BaNiF4 19 ; in Ref. 19 it was 
found that at U = 5 eV, the canting varies from 2° to 3° 
when J is varied from eV to 1 eV. In both LiNiP04 and 
BaNiF4 the Ni ion is divalent, with a d 8 configuration, 
and octahedrally coordinated. To investigate the gener- 
ality of this behavior, we next consider the case of the 
canted-spin antiferromagnet NiF 2 , in which the Ni ion is 
in the same coordination environnement as in BaNiF4. 
Experimentally, NiF 2 has the spins aligned preferentially 
in the plane perpendicular to the c axis with a slight cant- 
ing from antiparallel alignment by an estimated ~0.5° at 
low temperatures 13 . Performing LSDA+Z7 calculations 
at the experimental volume and with U = 5 eV and J = 
eV we indeed obtain the easy axis perpendicular to the c 



axis and a small canting of 0.3°, in excellent agreement 
with the experiments. In contrast to the case of LiNiP04, 
however, we find that the amplitude of the canting angle 
is almost insensitive to the value of J with just a small 
tendency to be reduced when J increased. This insen- 
sitivity of the canting angle to the value of J in NiF 2 
can be understood from the fact that in this compound 
the magnetism is almost collinear, and therefore the off- 
diagonal elements of the occupation matrix, and n^, 
are close to zero. Inspection of Eqs. 3 then shows that 
the effect of J is reduced largely to the diagonal part of 
the potential where the U parameter is dominant. 

To summarize our findings for the Ni-based com- 
pounds, in cases where the experimental canting is large 
(2-3°) we find a strong J-dependence of the canting an- 
gle, which increases with increasing J; when the canting 
is weak experimentally the J-dependence is much weaker. 




FIG. 2. Magnetocrystaline anisotropy energy versus the J parameter of (a) FeF2 (Experimental value from Ref. 20), (b) N1F2, 
(c) MnF2 (Experimental value from Ref.21) and (d) C0F2 ("sc" are calculations with Co semi-cores while "no sc" are calculations 
without Co semi-cores). The MCAE reported here is the energy between the a and c orientation of the spins, the energy of the 
c orientation is taken as reference. 



Next we analyse the effect of J on the behavior on 
the corresponding divalent iron compounds. We begin 
with LiFeP04, which is known experimentally to be a C- 
type AFM with an easy axis along the b direction and no 
observed canting of the spins 22 ' 23 (C y ground state with 
mmm' magnetic point group). Our calculations within 
the LSDA+/7 functional at the commonly used values of 
U = 4 eV and J = eV for Fe 2+ yield the correct C- 
type AFM order but find the easy axis incorrectly along 
the a direction. Now we switch to J 7^ eV and report 
in Fig. l.c the MCAE between the b and a directions, 
calculated by turning all the spins homogenously from 
the C y to the C x direction. We find that the MCAE is 
approximately linear with J, but with rather dramatic 
qualitative dependence: while at J = eV the easy axis 
is along the a direction (negative MCAE) the MCAE is 
almost reduced to zero around J = 0.5 eV and the easy 
axis changes to the b direction for J > 0.5 eV (positive 
MAE). To reproduce the experimental easy axis (C y ) a 
value of J greater than 0.58 eV is required. In the cases 



where the correct easy axis is reproduced (C y ) we do not 
observe any canting of the spins, in agreement with the 
experimental magnetic point group mmm\ 

As a second example with Fe 2+ , we analyse the effect of 
J on the MCAE of FeF 2 . Experimentally FeF 2 is known 
to have its spin magnetization parallel to the tetragonal c 
axis with a rather large MCAE of about +4800 /ieV 24 ' 25 . 
In Fig. 2. a we report the LSDA+/7 MCAE energies with 
respect to J at four different values of U (3, 4, 5 and 6 
eV). All the calculations with J = eV give the wrong 
easy axis (spins are perpendicular to c) with a huge error 
in the MCA energy (MCAE from -16000 to -26000 /ieV 
for U going from 3 to 6 eV). Increasing the value of J in 
the range of 0-0.5 eV has the tendency to strongly reduce 
this error with a linear increase of the MCAE with J as 
we found above for LiFePC^. However beyond J 0.5 
the increase of the MCAE is reduced and the evolution 
becomes more complex with the appearance of two max- 
ima before a drastic decrease beyond J c± 1.3 eV. The 
correct easy axis (MCAE> 0) is only obtained for a very 
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small range of U and J values, and the amplitude of the 
MCAE is correct over an even smaller range. This J de- 
pendence of the MCAE is again consistent with Eqs.3-6. 
From Eq.4 it is clear that when changing the orienta- 
tion of the spins from the z axis to the x or y axis the 
off-diagonal parts of Eq.4 become non-zero resulting in a 
direct effect of J on the MCAE from Eqs.6. 

We also performed the same analysis of the MCAE 
for NiF 2 (Fig.2.b), MnF 2 (Fig.2.c) and CoF 2 (Fig.2.d). 
MnF 2 and CoF 2 have the same easy axis as FeF 2 while 
NiF 2 has its easy axis perpendicular to the c direction. 
The easy axis is well reproduced for all three compounds 
at J = eV. As for FeF 2 , the amplitudes of the MCAE 
depend strongly on J but with a completely different 
trend in each compound. For MnF 2 and FeF 2 the experi- 
mental value can be reproduced by adjusting the values of 
U and J. In the case of CoF 2 and NiF 2 no experimental 
values are available. For CoF 2 we also performed calcula- 
tions with and without Co semi-cores states (Fig.2.d) and 
find a strong difference in the magnitude of the MCAE for 
the two cases. For FeF 2 we also performed calculations 
within the GGA functional (black pentagons in Fig. 2. a) 
and obtained a completely different J dependence than 
those calculated with the LDA functional. These com- 
parisons illustrate the difficulty of extracting a general 
rule about the J dependence of the MCAE. 

Our results reveal a problem with the predictability of 
the LSDA+Z7 method for non-collinear magnetic materi- 



als: A strong dependence of the MCAE and spin canting 
angles on the values of U and particularly J that are 
used in the calculation. Since properties such as magne- 
tostriction, piezomagnetic response, magnetoelectric re- 
sponse and exchange bias coupling are directly related to 
MCAEs and spin canting, it is of primary importance 
to reproduce these quantities accurately. At the mo- 
ment, the most reliable, although not entirely satisfac- 
tory, option appears to be a fine tuning of the U and J 
parameters by adjustment to reproduce experimentally 
measured anisotropics and canting angles; there is some 
evidence to suggest that properties such as magnetoelec- 
tric responses are then in turn well reproduced 26 . Future 
studies might explore methodologies for self-consistent 
calculation of the J parameter, or the predictions of new 
descriptions of the exchange and correlation such as the 
hybrid functionals 27 . On the flip side, it is clear that non- 
collinear magnetic systems provide a challenging case for 
testing the correctness of new exchange correlation func- 
tionals within the density functional formalism. 
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